Strange and charm quark-pair production 
in strong non-Abelian field 



Peter Levai 

KFKI RMKI Research Institute for Particle and Nuclear Physics, 
P.O. Box 49, Budapest 1525, Hungary 

E-mail: plevai@rmki.kfki.hu 

Vladimir V. Skokov 

Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 
Dubna, 141980, Russia 

GSI, Planckstrafie 1, D-64291 Darmstadt, Germany 
E-mail: V.Skokov@gsi.de 

Abstract. 

We have investigated strange and charm quark-pair production in the early stage 
of heavy ion collisions. Our kinetic model is based on a Wigner function method 
for fermion-pair production in strong non-Abelian fields. To describe the overlap of 
two colliding heavy ions we have applied the time-dependent color field with a pulse- 
like shape. The calculations have been performed in an SU(2)-color model with finite 
current quark masses. For strange quark-pair production the obtained results are close 
to the Schwingcr limit, as we expected. For charm quark the large inverse temporal 
width of the field pulse, instead of the large charm quark mass, determines the efficiency 
of the quark-pair production. Thus we do not observe the expected suppression of 
charm quark-pair production connecting to the usual Schwinger-formalism, but our 
calculation results in a relatively large charm quark yield. This effect appears in 
Abelian models as well, demonstrating that particle-pair production for fast varying 
non-Abelian gluon field strongly deviates from the Schwinger limit for charm quark. 
We display our results on number densities for light, strange, charm quark-pairs, and 
different suppression factors as the function of characteristic time of acting chromo- 
electric field. 



PACS numbers: 24.85. +p,25.75.-q, 12.38.Mh 
1. Introduction 

In the transport models theoretical descriptions of particle production in high energy 
pp collisions are based on the introduction of chromoelectric flux tube ('string') models, 
where these tubes are connecting quark and diquark constituents of colliding protons [1]. 
However, at RHIC and LHC energies the string density is expected to be so large 
that a strong collective gluon field will be formed in the whole available transverse 
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area. Furthermore, the gluon number will be so high that a classical gluon field 
as the expectation value of the quantum field can be considered in the reaction 
volume [2| We have investigated quark-pair production and determined particle 
spectra in time-dependent external U(l) and SU(2) chromo-electric fields |11[5]. In this 
paper, we describe strange and charm quark-pair production, and make a calculations 
of corresponding suppression factors for SU(2) gauge field. The results of solving similar 
problem for U(l) gauge field can be found in |6j. 

An alternative approach, that takes into account space inhomogeneities, was 
considered in (TIE]. However, it is worthwhile to mention that in contrast to the main 
idea of [3 [H] , where pairs production is directly calculated by numerical integration of 
a Dirac equation, our approach based on solving a kinetic equation for an "observable" 
Wigner function (or, finally, distribution function) providing to a considerable extent an 
intuitive insight to the physical problem. The next advantage of the current approach 
that it is not so highly computer demanded, thus allows to obtain detailed information 
about created particles. 



2. The kinetic equation for the Wigner function 



F,.,{iy,7V} 



The equation of motion for color Wigner function in gradient approximation reads flUi [9] : 

W + |^(4WFo.} + 2{F,.,[iy,7V]}- 
= zfc,{7°7\ W} - zm[70, W] + ig [A, , [7°7\ W] 

The color decomposition of the Wigner function with SU(A''c) generators in 
fundamental representation is given by 



1) 



w = w' + we 



12 A^^ 



1 



(2) 



where is the color singlet part and W"" is the color multiplet components. It is also 
convenient to perform spinor decomposition separating scalar a, vector 6^, tensor c^i,, 
axial vector d^ and pseudo-scalar parts e: 

The asymmetric tensor components of the Wigner function can be decompose into axial 



and polar vectors cj = c^^ and C2 



Cuip correspondingly. 



3. Kinetic equation with SU(2) color isotropic external field 

After decomposition the equations for the Wigner function in the case of pure 
longitudinal external SU(2) color field. A" = (0, 0, 0, with fixed color direction 

A1 = A^n"-, where ri^n" = 3 and dtu"' = [5J, we obtain the following system of 
equations for singlet components 

dta' + ^Et^a^ = -Akcl, (4) 
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d,e^ + ^K^e'' = -Ak4 - 2mdl (5) 

^*^o + T^^^I:^-0, (6) 

dth^ + = 2[fc X d^] + 4mct, (7) 

d,dl + ^i^t^dl = 2me^ (8) 

d,d^ + ^Et^d^ = 2[kxb% (9) 

5,c^i + = ci'k - m^,^ (10) 

dtC^2 + ^E:^4 = e'k; (11) 
and multiplet components 

dta^ + gE:-^a' = -Akct, (12) 

9ie<^ + gEt^e' = -4k4 - 2mdl (13) 

d,hl + gE:^¥, = Q, (14) 

5*6^ + gEt^b' = 2[k X d^] + Amcl, (15) 

dtdt + gEt^d',5'' = 2me'', (16) 

9,d^ + (7E:^d^ = 2[fcx6-], (17) 

dtcl + gEl^cl = a'^k - mb\ (18) 

d,cl + gEl^cf, = e''k, (19) 

where SU(2) triplet components of the Wigner function are defined by a"^ = a^n"". 



The distribution function of quarks (antiquarks) is defined by the components 
a',b' 0: 

^ /, N ^/ , N ma^ikA) + kb^ikA) 1 , /—^ ^ 

f,{k,t) = fq{~k,t) = -^^^ ^ + 2' ^ " Vfc2 + m2.(20) 

Thus to obtain the quark distribution function, the scalar a, vector 6^, axial vector d^, 
and axial tensor Ci components of the Wigner function are required, only. 
The initial conditions for the Wigner function in vacuum reads 0: 
1 m ^ 1 fc 



Strange and charm quark-pair production in strong non-Abelian field 



4 



and zero initial conditions for the rest components of Wigner function. Considering 
symmetry of initial condition and performing the vector decomposition, 



v = v^n + v±- \-v^[nx—], (22) 

we obtain the following equations for singlet components (c = ci to simplify reading): 

dta' + = -4(fc.c: + fc^cl), (23) 

+ = ^^^^^ + ^"""^ ^^^^ 

dtbl + ^K^K = -2Kdl + 4mcl, (25) 

+ = ^^^^^^ " ^^^^ 

5*c: + ^i?:^c: = a^A;,-m6:, (27) 

dtcl + ^Et^cl = a'k^ - mbl ; (28) 
and for multiplet components: 

dta' + gEt-^a' = -A{k,cl + fc^cl), (29) 

dtK + = 2fc±< + 4mc:, (30) 

5*61 + gEt^h^ = -2k,dl + Amcl, (31) 

dtdl + gEt^dl = 2{k,bl - k^bt), (32) 

dtcl + gEl^cl = a'^k. - mbt, (33) 

dtcl + gEl-^cl = a'^k^ - mbl . (34) 



The axial part of vector bj. and tensor components 0^, longitudinal d^ and transverse d± 
parts of axial vector components do not contribute to the evolution of the distribution 
function and thus are not considered. 



4. Numerical results and discussions 



In Ref. [5] we have solved the above equations and described the time evolution of the 
quark distribution functions to obtain the longitudinal and transverse quark spectra. In 
contrast to Ref. [5], in the current paper we focus on the integrated particle yields. 

In the numerical calculation we have used the following parameters: the maximal 
string tension Eq = 1 GeV/fm; coupling constant g = 2; the current quark masses niu^d = 
8 MeV, rUs = 150 MeV, = 1200 MeV for light, strange and charm quark, respectively. 
The particle production is ignited by a pulse-like field E^{t) = Eq ■ 1 — tanh^(t/r 
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Figure 1. Left panel: the total quark-pair number densities for different flavours, 

1/2 

ni{t), as a function of time t for short pulse width tEq = 0.1 (solid lines) and long 

1 /2 

pulse width tEq — 0.5 (dashed lines). Right panel: the total quark-pair number 
densities at the final state, n(t 3> t), as a function of pulse width r. 



which is characterized by the amphtude of the pulse Eq and its temporal width r. In 
this treatment particles are produced and absorbed by the field pairwise. 

The ratio of number densities of heavy quark-pairs, e.g. strange, to light quark- 
pairs (u, d-quarks) is widely known as a suppression factor. In our model it is defined 
in the asymptotic future (c.f. [Hj), t ^ r, as 

7 = ^lim riheavy (t) I niight (t ) , (35) 

where Uqit) is the number density of corresponding quark-pairs given by 

d^k 



ng{t) 



with degeneracy factor u = 2 (spin) x 2 (quark — antiquark) x A^^ (color). 



(36) 
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Figure 2. Left panel: the pulse width, r, dependence of the suppression factor 7. 
Arrow indicates the Schwinger limit for strangeness suppression factor. Right panel: 
the quark mass, m, dependence of the suppression factor 7 at different pulse width. 
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In Fig. [T] the time evolution of quark-pair number densities, displayed 

1 /2 

for different pulse widths, tEq = 0.1 and 0.5. For short pulse width the quark-pair 
number densities are comparable with each other (solid lines). In this case the particle 
production happens during the whole evolution of the field. In contrast to this, for long 
pulse, the number densities of produced charm quark-pairs becomes negligible in the 
final state, because charm-pair production is balanced by subsequent absorption by the 
field. This dependence on the pulse width is also demonstrated on the right panel of 
Fig. [H This figure clearly displays that charm quark-pair production is substantially 
enhanced in the cases of short pulse widths, and this enhancement has a maximum 
at r ~ O.Ia/Eq. In opposite to the heavy charm quark, light and strange quark-pair 
productions are increasing with the pulse width, without any local maximum. 

We have investigated the suppression factor and its dependence on pulse widths and 
quark masses. Fig. [2] summarizes our results. On the left panel the dependence on the 
pulse width is displayed. The strange to light ratio has a weak dependence on the pulse 
width, its value is approaching slowly the asymptotic value of Schwinger limit (0.84) 
from below, similarly to U(l) gauge field p|. For charm quark this Schwinger limit is 
negligibly small, which value is reproduced by our numerical calculation for very long 
pulse width. On the other hand, at short pulse widths, the relative charm production 
is surprisingly large, which does not follow any earlier expectation. Considering charm 
to light and charm to strange ratios, only a slight difference can be seen between them. 
On the right panel we display the quark mass dependence of the suppression factor for 
different pulse widths. For short pulse width the suppression factor is decreasing almost 
linearly with increasing quark mass value. For large pulse widths we can see a very fast 
(~ exp{— m^/E'o}) drop, which is consistent with the Schwinger formula. 

5. Conclusion 

We have calculated light, strange and charm quark-pair production in time-dependent 
SU(2) non-Abelian field. Applying a pulse-like time evolution and investigating the 
influence of pulse width, we observed that hght and strange quark-pairs are produced 
as we expected, approaching the Schwinger limit. Charm quark-pairs followed this 
behaviour for large pulse widths. However, for short pulses we did not see the expected 
charm suppression, connected to the large charm quark mass. Indeed, the large value 
of inverse temporal width of the pulse, overwhelming the mass of the heavy quark, 
1/r ^ nic, determines the quark-pair production. This finding could indicate the 
formation of collective gluon field via enhanced heavy quark-pair production at RHIC. 
The issue of quantitative calculation of particle suppression factors for different quarks 
and comparison with the existing models will be addressed elsewhere. 
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